Lifetimes of Shockley electrons and holes at the Cu(lll) surface 
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A theoretical many-body analysis is presented of the electron-electron inelastic lifetimes of Shock- 
ley electrons and holes at the (111) surface of Cu. For a description of the decay of Shockley states 
both below and above the Fermi level, single-particle wave functions have been obtained by solving 
the Schrodinger equation with the use of an approximate one- dimensional pseudopotential fitted to 
reproduce the correct bulk energy bands and surface-state dispersion. A comparison with previous 
calculations and experiment indicates that inelastic lifetimes are very sensitive to the actual shape 
of the surface-state single-particle orbitals beyond the F (kn = 0) point, which controls the coupling 
between the Shockley electrons and holes. 

PACS numbers: 71.10.Ca,72.15.Lh,73.20.At,68.37.Ef 
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I. INTRODUCTION 

A variety of metal surfaces, such as the (111) surfaces 
of the noble metals, are known to support a partially oc- 
cupied band of Shockley surface states with energies near 
the Fermi level^ whose dynamics have been the subject 
of long-standing interesti 2 * 3 ^ In particular, the lifetimes 
of excited holes at the band edge (k|| = 0) of these sur- 
face states have been investigated with high resolution 
angle-resolved photoemission (ARP^SiS and with the 
use of the scanning tunneling microscope (STM). 9 STM 
techniques have also allowed the determination of the 
lifetimes of excited Shockley and image electrons over a 
range of energies above the Fermi level>i2iii 

Many-body calculations of the electron-electron (e-e) 
inelastic lifetimes of excited holes at the surface-state 
band edge of the (111) surfaces of the noble metals 
Cu, Ag, and Au, which were based upon the G°W 
one-dimensional scheme that had been introduced by 
Chulkov et al. 12 to describe the lifetimes of image- 
potential statesji 3 " showed considerable agreement with 
experiments^ These calculations were then extended 
to treat the case of excited surface-state and surface- 
resonance electrons above the Fermi levelSii In order to 
account approximately for the potential variation in the 
plane of the surface, the original one-dimensional model 
potential, which had been introduced to describe surface 
states at the T point, was modified along with the in- 
troduction of a realistic effective mass for the dispersion 
curve of both bulk and surface states. Within this model, 
however, all Shockley states have the same effective mass, 
so the projected band structure is not correct, especially 
at energies above the Fermi level, as shown in Fig. ^ 

In this paper, we present a new approach that although 
at the r point is less sophisticated than the model used 
in Refs. l3T5lll^lr^ (i) basically reproduces the surface- 
state probability density of Refs. Il4ll5ll6ll7l at the band 
edge of the surface-state band in Cu(lll) and (ii) has 
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FIG. 1: Dispersion of the Cu(lll) Shockley surface state 
(thick solid line), as obtained from three-dimensional ab initio 
calculations. Shaded areas represent areas outside the band 
gap, where bulk states exist. The thin solid and dotted lines 
represent approximate energy dispersions of the Shockley sur- 
face state and the bottom of the projected band gap, respec- 
tively, as obtained from e + k 2 /2m with e = — 0.44eV and 
m = 0.42 in the case of the surface state (thin solid line) and 
e = — 1.09eV and m = 0.22 for the bottom of the projected 
band gap (thin dotted line). 



the merit that it reproduces, through the introduction 
of a k|| -dependent one-dimensional potential, the actual 
bulk energy bands and surface-state energy dispersion 
of Fig. 1, thereby allowing for a realistic description of 
the electronic orbitals beyond the T point. Adding the 
contribution from electron-phonon coupling, 1 — which is 
particularly important at the smallest excitation ener- 
gies, our calculations of the lifetime broadening of excited 
Shockley electrons and holes in Cu(lll) indicate that (i) 
there is good agreement with experiment at the surface- 
state band edge and (ii) at energies above the Fermi level 
the lifetime broadening is closer to experiment and very 



sensitive to the actual shape of the surface-state single- 
particle orbitals beyond the T point. 

Let us consider a semi-infinite many-electron system 
that is translationally invariant in the plane of the surface 
(normal to the z axis) . The decay rate or inverse lifetime 
of a quasiparticlc (electron or hole) that has been added 
to the system in the single-particle state e lki ' r ip kiyEi (z) 
of energy Ei is obtained as follows (we use atomic units, 
i.e., e 2 — h = m P = 1)^ 



^ki,Ei = T'2JdzJ dz'ifjl i Et (z)lmJ:(z,z';k i ,E i )i} k ^E z {z / ), 

(1) 

where the =F sign in front of the integral should be taken 
to be minus or plus depending on whether the quasi- 
particle is an electron (Ei > Ep) or a hole (Ei < Ep), 
respectively, Ep is the Fermi energy, r and represent 
the position and wave vectors in the plane of the surface, 
and z'\ h, Ei) is the nonlocal self-energy operator. 

To lowest order in a series-expansion of the self-energy 
E(z, z'\ hi, Ei) in terms of the energy-dependent screened 
interaction W(z, z'\ k, E^jSLSi and replacing the interact- 
ing Green function G(z, z'; fc, E) by its noninteracting 
counterpart, one finds the following expression for the 
imaginary part of the so-called G°W self-energy: 



lmS(z,z';ki,Ei) = — / dE 

7T JO 

xImG°(z, z 1 ; k f ,E f ) JmW(z, z'\ q, E). 



dq 



(2) 



Here, kf — ki — q, Ef = Ei — E, and G°(z, z'; k, E) is 
the noninteracting Green function 



G°(z,z';k,E) 



^k,E^k,E ( Z ) 



(3) 



where z < (z > ) is the lesser (greater) of z and z 1 , and 
[f,g](z) is the wronskian 



[f,g](z)=f(z)g'(z)-f'(z)g(z). 



(4) 



The functions ip^ E ( z ) are solutions of the single-particle 
Schrodinger equation 

- (1/2) + V k (z) ^ E {z) = {E- k 2 /2) ^ k<E (z) (5) 

regular at ±oo, with V k (z) being a momentum-dependent 
one-dimensional effective potential that we fit to the pro- 
jected surface band structure. We use 



V k (z) 



U k + 2V k cos(2icz/a s ), z < z k 



(6) 



z > z k , 



where U k and V k are fitted to the bulk energy bands 
(which we have obtained from three-dimensional ab ini- 
tio calculations), a s = 2.08 A represents the interlayer 
spacing, $ = 4.94 eV is the experimentally determined 
work function, and the matching plane z k is chosen to 
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FIG. 2: Probability-density \if) k ,E\ 2 of the Shockley surface 
state at the center of the surface Brillouin zone of Cu(lll), 
as obtained within two different models: the model pre sented 
here (dashed line) and the model reported in Ref. [l^ (solid 
line). The dotted line represents the probability-density that 
we obtain at k = 0.2 a.u. Full circles represent the atomic po- 
sitions of Cu in the (111) direction. The geometrical (jellium) 
electronic edge (z = 0) has been chosen to be located half an 
interlayer spacing beyond the last atomic layer. 



give the correct surface-state dispersion represented in 
Fig. 1 by a thick solid line^ 

In the random-phase approximation (RPA)^ the 
screened interaction W(z, z';q, E) is obtained from the 
knowledge of the noninteracting density-response func- 
tion x°( z j z '] E) by solving the following integral equa- 
tion 

W(z, z'\ q, E) — v(z, z'\ q) + J dz\ J dz 2 
xv{z, zi;q) x°(zi, z 2 ; q, E) W(z 2 , z'\ q, E), (7) 

where v(z, z'; q) represents the two-dimensional Fourier 
transform of the bare Coulomb interaction. The re- 
sults presented below have been obtained by using in 
X (z, z';q,E) the eigenfunctions and eige nvalues of the 
one-dimensional Hamiltonian of Ref. [l2j with all effec- 
tive masses set equal to the free-electron mass. We have 
also used the eigenfunctions and eigenvalues of a single- 
particle jellium-surface Kohn-Sham Hamiltonian (in the 
local-density approximation) with r s = 2.67^* and we 
have found that the surface-state lifetimes are rather in- 
sensitive to whether one or the other choice is employed. 

The abrupt step model potential of Eq. @ , which does 
not account for the image tail outside the surface, could 
not possibly be used to describe image states. However, 
Shockley surface states are known to be rather insensi- 
tive to the actual shape of the potential outside the sur- 
face; indeed, the model potential of Eq. is found to 
yield a surface-state probability-density \ip ki ,Ei\ 2 at the 
band edge of the Shockley surface-state band of Cu(lll) 
(ki = 0) that is in reasonably good agreement with the 
more realistic surface-state probability density used in 
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TABLE I: G°W decay rates, in linewidth units (meV), of an 
excited hole at the band edge of the Shockley surface-state 
band of Cu(lll) (Ei = -0.44 eV and h = 0). r intcr and r lntra 
represent interband and intraband contributions to the e-e de- 
cay rate r e - e - The total decay rate Ttotai includes the e-ph 
decay rate of 7 meV reported in Ref. ITsl The full calculation 
represents the result we have obtained by using in Eq. the 
actual fc-dependent model single particle wave function ipk,E- 
The approximate calculation represents the result we have ob- 
tained by replacing all surface-state wave functions with k / 
by that at the T point . The third row represents the calcula- 
tions reported in Refs. ll4jl5l The experimental linewidth has 
been taken from the STM measurements reported in Ref. Hil 



Calculation 


Tintcr 


Tintra 




Ttotal 


Full 


10 


9 


19 


26 


Approximate 


10 


23 


33 


40 


Refs. 14.15 


6 


19 


25 


32 


Experiment 








24 



Refs. IllllilllllO as shown in Fig. 2. Both probability- 
densities coincide within the bulk, although our approx- 
imate probability-density appears to be slightly more 
localized near the surface, as expected. 2 - Nevertheless, 
we find that decay rates of an excited hole at f based 
on the use of these two models to describe the wave 
function ipki, Ei entering Eq. Q agree within less than 
1 meV. Differences between our new calculations, which 
are based on the use of the fc-dependent model potential 
of Eq. |JB}, and those reported previously ; 14 ! 15 ! 16 ! 17 are 
due to our more realistic description of the band struc- 
ture and surface-state wave functions beyond the T point. 

First of all, we consider an excited hole at the band 
edge of the Shockley surface-state band of Cu(lll), i.e., 
with Ei = -0.44 eV and k t = (see Fig. 1). The decay 
of this excited quasiparticle may proceed either through 
the coupling with bulk states (interband contribution) 
or through the coupling, within the surface-state band 
itself, with surface states of different wave vector k par- 
allel to the surface (intraband contribution). In order to 
investigate the impact of the actual shape of the surface- 
state wave functions with k ^ on the decay of the 
surface-state hole at T, we have compared in Table I 
the decay rates that we have calculated either by using 
in Eq. (|3J) the actual fc-dependent model wave function 
V-'fc.-E (full calculation) or by replacing all surface-state 
wave functions ipk.E with k =/= by that at the T point 
(approximate calculation). This comparison shows that 
the penetration of the actual fc-dependent surface-state 
wave functions ipk,E being larger than at T (compare the 
dashed and dotted lines of Fig. ^) yields a reduction in 
the decay rate from 33 meV to 19 meV, which is due to 
the fact that the coupling of the surface-state hole at f 
with actual surface states of different wave vector k (in- 
traband contribution) is smaller than the coupling that 
would take place with surface-state orbitals that do not 
change with k . The difference between our predicted 
surface-state lifetime broadening of 19meV and that re- 



ported before (t _1 = 25meV)^i^ is entirely due to our 
more accurate description of (i) the projected band struc- 
ture and (ii) the wave-vector dependence of the surface- 
state wave functions ipk,E(z) entering the evaluation of 
the Green function of Eq. J3J). 

We have also carried out a full calculation of the decay 
of an excited hole at T but replacing the actual surface- 
state wave vector kf entering Eq. (0) by the wave vec- 
tor that would correspond to a parabolic surface-state 
dispersion of the form dictated by the thin solid line of 
Fig-ffl an d we have found that the linewidth is reduced (as 
expected, since the parabolic dispersion results in fewer 
final states) by less than ~ 1 meV. However, if one fur- 
ther replaces our wave- vector dependent surface-state or- 
bitals entering Eq. J2J) by their less accurate counterparts 
used previouslyji*i^ the lifetime broadening is increased 
considerably (from 19 to 25meV), showing the impor- 
tant role that the actual coupling between initial and 
final states plays in the surface-state decay mechanism. 

Our model, which correctly reproduces the behaviour 
of s-p valence states, does not account for the presence 
of d electrons with energies a few electronvolts below the 
Fermi level. The screening of d electrons is known to play 
a crucial role in the decay mechanism of bulk states^ 
However, in the case of Shockley holes, whose decay is 
dominated by intraband transitions that are associated 
with very small values of the momentum transfer, the 
screening of d electrons is expected to reduce the life- 
time broadening only very slightly^ and will not be in- 
cluded in the present work. Adding to our estimated e-e 
linewidth of 19 meV the elect ron-phonon (e-ph) linewidth 
of 7eV reported in Ref. fl8L we find r to tai = 26meV 
in close agreement with the experimentally measured 
linewidth of 24 meV, as shown in Table I. 

In Fig. 3, we show our full calculation (thick solid line) 
of the inelastic linewidth (r to tai = Tc-c + To-ph) of ex- 
cited Shockley holes and electrons in Cu(lll) with en- 
ergies Ei below and above the Fermi level. Also shown 
are separate interband and intraband contributions to 
the linewidth (inset), the approximate calculation that 
we have carried out by replacing all surface-state wave 
functions with k ^ by the surface-state wave function 
at r (thin solid line), and the calculations reported in 
Ref. [l6| (thin dashed line). The lifetime broadening of 
excited Shockley electrons in Cu(lll) was studied with 
the STM by Biirgi et al.ML However, it has been argued 
recently^ that due to an error in the identification of the 
phase coherence length in the measured quantum inter- 
ference patterns the linewidths reported in Ref. should 
be doubled. These doubled values have been represented 
in Fig. 3 by solid squares. The experimentally deter- 
mined inelastic linewidth of an excited Shockley hole at 
the F points is represented by a solid circle. 

Fig. 3 shows that a correct description of the wave- 
vector dependent surface-state wave functions reduces 
the coupling of holes and electrons within the Shock- 
ley band, thereby bringing the lifetime broadening into 
closer agreement with experiment. A comparison be- 
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FIG. 3: Scaled total linewidth T totsl /(Ei - E F ) 2 of Shockley 
surface-state electrons and holes in Cu(lll) (including the 
e-ph contribution r e _ p h = 7meV) versus the surface-state 
energy Ei, as obtained from Eqs. 0-0 (thick solid line), 
from Eqs. (0-0 but replacing all surface-state wave func- 
tions with k 7 ^ by those at the f point (thin solid line), and 
from Ref. [ly (thin dashed line) . The dotted line represents 
the energy-independent e-ph contribution r e _ p h = 7meV. 
The dotted line with squares represents the STM measure- 
ments reported in Ref. YWL multiplied by a factor of two to 
correct for an error in the phase coherence length used in 
that work. The solid circle represents the experimentally 
determined linewidth of an excited hole at the f point, as 
reported in Ref. Q. Separate unsealed interband (Anter) 
and intraband (rj ntr a) contributions to the total linewidth 
(r to tai = Tintra + Tintor + IV ph ) are represented in the inset 
by dashed and dottel lines, respectively. 

tween our full calculations and experiment shows that 
there is close agreement at the surface-state band edge 



(at E — Ep = — 0.44eV) and there is also reasonable 
agreement at energies above the Fermi level. At ener- 
gies where the surface-state band merges into the contin- 
uum of bulk states, however, our calculated linewidths 
are still too low, which might be a signature of the 
need of a fully three-dimensional description of the sur- 
face band structure. We also note that differences be- 
tween the calculations reported here and those reported 
previousljii 6 indicate that inelastic lifetimes are very sen- 
sitive to the actual shape of the surface-state single- 
particle orbitals beyond the T point. The linewidths re- 
ported here are smaller (larger) for excited holes (elec- 
trons) below (above) the Fermi level, thus bringing the 
theretical predictions closer to experiment. 

In summary, we have presented a new G°W one- 
dimensional scheme to calculate the inelastic lifetime 
broadening of excited Shockley electrons and holes in 
Cu(lll), which is based on a realistic description of the 
projected bulk energy bands and the surface-state or- 
bitals beyond the T point. Adding the contribution from 
electron-phonon coupling^ which is particularly impor- 
tant at the smallest excitation energies, our calculations 
indicate that there is reasonable agreement with experi- 
ment, especially at low excitation energies. The screening 
of d electrons, not included in this work, is expected to 
reduce the lifetime broadening only very slightly, at least 
at the hot-electron energies nearest to the Fermi level. 
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